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Abstract 

We present a search for a pentaquark decaying strongly to H~7r~ in 7^ collisions 
at a center-of-mass energy up to 25 GeV/c 2 . Finding no evidence for such a state 
in the mass range of 1480 MeV/c 2 to 2400 MeV/c 2 , we set limits on the yield and 
on the cross section times branching ratio relative to H*(1530)°. 
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1 Introduction 

The existence of bound multiquark states like QQqq and the H dihyperon 
were first proposed by Jaffe [1,2,3] in 1977. Then 20 years later Diakonov et 
al. [4] proposed the existence of four quarks and one antiquark confined in 
a low-mass anti-decuplet configuration. In their calculations Diakonov made 
several predictions of masses and widths of exotic baryonic states such as the 
mass of the lightest state B + (previously called Z + ) at about 1530 MeV/c 2 . 
There were also predictions of decay modes: + — ► pK° s , 6 + — > nK + and 
ET - H-7T- 0. 

1 See http://www-focus.fnal.gov/authors.html for additional author information. 

2 The is also known as the 0(1860) . 
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On the experimental side, the year 2003 was the beginning of "pentaquark 
observations." The B + at about the predicted mass was the first candidate 
as noted in the PDG2004 [5], but searches by higher statistics experiments 
yielded negative results as noted by R. Schumacher [6]. 

In the case of the doubly strange pentaquark the only positive evidence is 
from NA49 [7]. Negative results were obtained by every other search: HERA-B 
[8], ALEPH [9], WA89 [10], HERMES [11], BABAR [12], ZEUS [13], COM- 
PASS [14], E690 [15], as well as preliminary results from CDF [16]. 

This letter describes a search for the — > 5~7r~ pentaquark candidate^! in 
electromagnetic interactions and extends the search for the singly strange + 
discussed in [17] to the doubly strange state considered here. 



2 Event reconstruction and selection 

The FOCUS experiment took data during the 1996-7 fixed-target run at Fer- 
milab. A photon beam obtained from bremsstrahlung of 300 GeV electrons 
and positrons impinged on a set of BeO targets. The first element in the 
spectrometer was the silicon strip detector array: four triplets of silicon strip 
planes used for track finding and vertexing. Each triplet was comprised of 
three closely spaced parallel planes with the strip directions rotated to pro- 
vide a means for correlating hits in the three planes of the triplet and thus 
track coordinates. For most of the run, two pairs of silicon strips were also 
interleaved with the target segments for more precise vertexing [18]. Charged 
particles were tracked and momentum analyzed as they passed through one or 
two dipole magnets and three to five sets of multiwire proportional chambers 
with four planes each (as shown in Fig. 1). Three multicell threshold Cerenkov 
counters, two electromagnetic calorimeters, and two muon detectors provided 
particle identification. A trigger which required, among other things, >25 GeV 
of hadronic energy passed 6 billion events for reconstruction. 

The data used for this analysis come from a subset of FOCUS data which 
contain cascade candidates (H~ — > A°7r~ and Q~ — > A°K~). The cascade 
decays used in this analysis are those which occur downstream of the sili- 
con detector and with the A daughter being fully reconstructed through the 
decay A — > pn~ . The candidate is a reconstructed silicon track with di- 
rection and position consistent with the intersection of a reconstructed A and 
a multiwire chamber track. The A invariant mass is required to be between 
1.10 GeV/c 2 and 1.13 GeV/c 2 . The higher momentum track is chosen to be 
the proton because the small phase space of the A — > piT~ constrains the 



Charged conjugate states are implied unless explicitly stated otherwise. 
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Fig. 1. A schematic drawing in the bend view of the spectrometer of a 
H°(1530) — > H _ 7r + decay, with the S~ decaying downstream the silicon strip detec- 
tor, as described in the text. Only the front part of the spectrometer is displayed. 

proton to carry most of the momentum for the decays observed in the for- 
ward FOCUS spectrometer. For this analysis we identify S~ tracks when the 
matched A candidate and the matched pion track have an invariant mass 
within 20 MeV/c 2 of the nominal S~ signal peak (shadow region of fig. 2). 
The total sample yields approximately 625 000 S - — > A°7r~ signal events. A 
detailed description of the reconstruction of cascades and vees can be found 
in Ref. [19]. 




Fig. 2. The invariant mass plot represent the total sample of silicon tracks matched 
to a reconstructed A and a multiwire chamber track (vr— ). The shadow area is the 
region used to select H™ tracks for this analysis. 

Each S~ track is combined with good quality charged tracks to make a vertex 
to search for S*(1530)° — > H~7r + and H5 ~ — > H~7r~. This two track vertex must 
be well defined with a x 2 probability greater than 1%. A multitrack production 
vertex is nucleated around this two track vertex and must be within 2a of the 
2~7r vertex. The production vertex must have a x 2 probability greater than 
1% and must be within the target material or outside by no more than 3a. In 
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both cases, a is the calculated uncertainty on the vertex location or separation. 



A particle identification algorithm has been developed which combines data 
from all of the Cerenkov counters which the track passes through. This al- 
gorithm [20] returns negative log-likelihood (times two) values Wi for a track 
and hypothesis i G {e, n, K,p} based on the light yields in the phototubes cov- 
ering the Cerenkov cone of the track. A Cerenkov cut, Wxcbi{e,K,p) ~ > —8, 
requires that the pion H~ daughter (tt~ track matched to 5~ track in Fig. 1) 
must not be strongly inconsistent with the pion hypothesis. For pions of the 
2~7T combination (n + track in Fig. 1) a Cerenkov cut, W m i n ( e! K, p ) — W n > —6, 
is used to reduce combinatorial background. 

Mass plots for S _ 7r + are shown in Fig. 3 and in Fig. 4 with an additional cut 
that the H _ 7r + momentum be greater than 25 GeV/c. The signal was best fit 
with a P-wave Breit-Wigner with an energy dependent width convoluted with 
a Gaussian for the detector resolution. The resolution, 2.92 MeV/c 2 , was ob- 
tained from a Monte Carlo simulation and the fitted width of the Breit-Wigner 
is shown in the figures, consistent with the widths of 8-10 MeV/c 2 quoted in 
the PDG[5]. The background was fit to the form aq b exp (cq + dq 2 + eq 3 + fq 4 ) 
where a — / are free parameters and q is the Q-value (invariant mass minus 
component masses). 
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Fig. 3. £*(1530)° fit with a P-wave Breit-Wigner and combinatorial background. 



3 Pentaquark search results 



The S~7r~ and S + 7r + invariant masses are plotted using the standard selec- 
tion criteria in Fig. 5. There are no significant differences between the two 
charge states so for the remainder of the analysis we combine the charge con- 
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Fig. 4. H*(1530)° fit with a P-wave Breit-Wigner for momentum greater than 25 
GeV/c and combinatorial background. 

jugate states. The combined sample with standard cuts and with an additional 
momentum cut of 25 GeV/c is plotted in Fig. 6. In this analysis we treat the 
two samples displayed in Fig. 6 separately, because the production mechanism 
utilized in Monte Carlo acceptance calculations is much better understood at 
the higher momenta. The same parameterization of the background is used in 
Fig. 6 as in Fig. 3 and 4. Such parameterization show acceptable fits to this 
smooth background as shown in Fig. 7. No evidence for a pentaquark near 
1860 MeV/c 2 or at any mass less than 2400 MeV/c 2 is observed. To set a limit 
on the yield we need to make some assumptions about the width of the state. 
We consider two cases: one with a natural width of and one with a natural 
width of 15 MeV/c 2 In the first case, the signal is fit with a Gaussian with a 
width given by the experimental resolution. In the second case, the signal is fit 
with a P-wave Breit-Wigner with an energy dependent width convoluted with 
the experimental resolution. The experimental resolution a is parametrized as 
a function of the invariant mass and we find a = —9.35 + 7.76m + 0.21m 2 
an adequate approximation, with a in MeV/c 2 and m is the mass in GeV/c 2 . 
With the additional momentum cut of 25 GeV/c the experimental resolution 
change to a = -7.99 + 6.29m + 0.62m 2 . 



A series of 921 fits to the observed S~7r _ mass plot were performed using 
the background and signal shapes described above for each assumed width. 
The signal mass is varied in 1 MeV/c steps from 1480 to 2400 MeV/c 2 and 
a binned log-likelihood fit using Minuit [21] is performed. The ±lcx errors 



NA49 detected this state with a width below the detector resolution of 18 
MeV/c 2 . Therefore the expected width of this state cannot exceed this value. A 
choice of natural width of 15 MeV/c 2 convoluted with FOCUS resolution in a range 
of 4-11 MeV/c 2 gives about 18 MeV/c 2 , the upper limit set by NA49. 
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Fig. 5. The invariant mass distribution of E - 7r - separated by charge (particle and 
antiparticle) . Standard cuts are applied (no momentum cut). The inset plot is the 
ratio of particle/antiparticle. 
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Fig. 6. The invariant mass distribution of 3~ir~ for particle and antiparticle com- 
bined. The solid line shows the result for the standard cuts and the dashed line is 
with the additional cut that the momentum is greater than 25 GeV/c. The inset 
plot is the ratio of both, note that the momentum cut affects primarily low mass. 
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Fig. 7. Fit to the invariant mass distribution of H~7r~ for particle and antiparticle 
combined. The left figure show the result for the standard cuts and the right figure 
with the additional cut that the momentum is greater than 25 GeV/c. 
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Fig. 8. Pentaquark yields and upper limits. Top (bottom) plots show results for 
a natural width of (15 MeV/c 2 ). The shaded region includes the la errors with 
the central value in the middle. The outer curves show the upper and lower 95% 
confidence limits. 

are denned as the point where Alog£ = 0.50 relative to the maximum log£, 
while continually adjusting the background parameters to maximize log£. 
The 95% CL lower limit is denned similarly with Alog£ = 1.92. Both are 
obtained using Minos [21]. The 95% CL upper limit is constructed as follows: 
The likelihood function £ versus yield is determined by maximizing log£ for 
many different (fixed) yields, allowing background parameters to float. The 
likelihood function is integrated from a yield of to oo to obtain the total 
likelihood. The 95% CL upper limit on the yield is defined as the point where 
95% of the total likelihood is between a yield of and the upper limit{£] The 
fitted yield, la errors, and 95% CL limits are shown in Figs. 8 and 9. Of the 
1842 fits, none of them finds a positive excursion greater than 5a. The only 
previous pentaquark observation was around 1860 MeV/c 2 . In this region, 
we find a small dip which is not statistically significant. The largest positive 
excursion occurs in the region where the background distribution peaks. 

To compare with other experiments, the limits on yield must be converted 
to limits on production times the (unknown) branching ratio. We choose to 

5 This definition of an upper limit is used rather than a counting based Feldman- 
Cousins type limit because errors are Gaussian for this large background. 



8 



I 
t 

1 

< 
w 
T 

E 




P>25GeV/c 



M(S 71 ) 



2.3 

GeV/c 



2.4 

2 



Fig. 9. Pentaquark yield and limits for standard cuts including the 25 GeV/c 
momentum cut. Top (bottom) plots show results for natural width of (15 MeV/c 2 ). 
The shaded region includes the la uncertainty with the central value in the middle. 
The outer curves show the upper and lower 95% confidence limits. 

normalize the production cross section to S*(1530)° because the recon- 



structed decay mode of the 5* (1530) 



S 7T 



is very similar, in terms of 



topology and energy release, to the state we are investigating. Thus, we at- 
tempt to determine 



Cr(»5 ~) ■ BR(~s- -> S-7T-) _ y(S 5 ~~) ■ gggT ~+ S-7T-) 6 S *(1 5 30)° 

o-(H*(1530)°) Y(S*(1530)°) £s , 

(1) 

All of the efficiencies include the reconstruction and selection efficiencies and 
corrections for unseen decays of parent particles. The 5* (1530)° efficiency in- 
cludes BR(S*(1530)° — > S _ 7T + ) = 0.66 and both efficiencies include the factors 
BR(H~ — > A°7r~) = 1 and BR(A — ► pn~) = 0.64. The two branching ratios 
which are common to both cancel. Determining reconstruction and selection 
efficiency (including acceptance) is described below. 

The FOCUS detector is a forward spectrometer and therefore acceptance de- 
pends on the momentum of the produced particle. The production characteris- 
tics of the pentaquark are the largest sources of systematic uncertainty in this 
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analysis. We choose a particular production model to obtain limits and pro- 
vide sufficient information about the experiment for other interested parties 
to obtain limits based on other production models. The production simula- 
tion begins with a library of e~ and e + tracks obtained from a TURTLE 
simulation [22] of the Wideband beam line. From this library, an individual 
track is drawn and bremsstrahlung photons are created by passage through a 
20% Xq lead radiator. Photons with energy above 15 GeV are passed to the 
Pythia [23] Monte Carlo simulation. The Pythia version we use is 6.127. 
The Pythia simulation is run using minimum bias events^] with varying 
energies!^] Options controlling parton distributions and gluon fragmentation 
were set to avoid heavy quark production^ Since Pythia does not produce 
pentaquarks, another particle must be chosen to represent the pentaquark. Ac- 
cording to the string fragmentation model, which is implemented in Pythia, 
the mass of the particle has the greatest effect on production and the number of 
quarks a particle has in common with the initially interacting hadrons is next 
in importance. The S*(1530)° particle is chosen to represent the production of 
a pentaquark. The S*°(ssm) can obtain at most 33% of the remaining quarks 
from the target nucleon valence quarks, while the (ddssu) can take 40%. 
The charge conjugate S*(1530)°(ssw) particles must obtain all quarks from 
the vacuum, while the ~(ddssu) can take 20% from the target nucleon. The 
mass of the particle chosen to re pre sent the pentaquark, 5* (1530)°, is set to 
the appropriate value in Pythia] 9 1 

To calculate the relative cross sections in Eq. 1 we need efficiencies for H*(1530)° 
2~7r + and — > 2~7r~. These efficiencies are obtained from the FOCUS 
Monte Carlo simulation. The dominant uncertainty in the efficiency determi- 
nation is the modeling of the production characteristics of the parent par- 
ticle. For the observed particle, S*(1530)°, we can compare the data and 
Monte Carlo directly and adjust the Monte Carlo simulation to produce the 
correct data distribution. Even this is not sufficient, however, because areas 
where the efficiency is zero cannot be accounted for. For S*(1530)°, we run a 
weighted Monte Carlo simulation which matches the Monte Carlo momentum 
distribution with the observed data momentum distribution in the region for 
which the acceptance is not zero. The dominant source of uncertainty for the 
S*(1530)° — > H~7T + efficiency is our lack of knowledge of the fraction of events 
completely outside of our acceptance (momentum less than 15 GeV/c). The 
weighted Pythia Monte Carlo predicts that 71% of the 5* (1530)° particles 
are produced with momentum less than 15 GeV/c. To obtain an estimate of 
the efficiency uncertainty, we assume that the number of particles with mo- 



6 Specifically MSEL=2 in the PYTHIA setup. 

7 Specifically MSTP(171)=1 in the Pythia setup. 

8 MSTP(58)=3 to produce light quarks (uds) only and MDME(156 — 160,1)=0 to 
limit gluon fragmentation into light quarks only in the Pythia setup. 

9 By setting PMAS(190, 1) in Pythia. 
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Fig. 10. Acceptance versus mass for pentaquark candidates. Upper (lower) curve is 
for a lo" pentaquark produced uncertainty. 

mentum less than 15 GeV/c can be off by as much as a factor of 2 (high or 
low). This leads to a relative uncertainty on the S*(1530)° efficiency of 45%. 
The — > 2~7T - efficiency is taken to match the uncertainty of a higher 
statistical mode + — > pK$ (~5%) [17] when using a substitute particle like 
H*(1530)° and E*(1385) + to generate pentaquarks with a momenta greater 
than 25 GeV/c. The S^ - — > H~7r~ efficiency versus mass (with no branch- 
ing ratio corrections) is shown in Fig. 10. The uncertainty in e~ — _>c-„- is 
approximately 10%. The relative uncertainty of the efficiency of an unknown 
particle (~10%) is indeed less than that for the high statistics normalizing 
modes (>40%) because the efficiency uncertainty of the high statistics modes 
reflects the lack of knowledge of production outside of our acceptance. It is 
reasonable to assume that discrepancies in the Monte Carlo simulation are 
similar for the signal mode and the normalizing mode since discrepancies are 
correlated and therefore adding the uncertainty to the signal mode would be 
double-counting. Note that the signal and normalizing efficiencies only appear 
ratio. 

We also report the relative cross sections in the region where our acceptance 
is good, that is for parent particle momenta greater than 25 GeV/c. This 
dramatically reduces the systematic uncertainties associated with the mea- 
surement. The uncertainty due to the production of S*(1530)° is minimal. 
The uncertainty in the H5 ~ efficiency is also reduced from approximately 10% 
to about 5% as shown in Fig. 11. The number of reconstructed S*(1530)° at 
momenta greater than 25 GeV/c is about 55 000, compared to a total sample 
of about 65 000 without the momentum cut. 

The upper limit on the yield was obtained by mathematically integrating the 
likelihood function from to infinity and then integrating from to 95% of 
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Fig. 11. Acceptance versus mass for pentaquark candidates. Lower (upper) curve is 
for la pentaquark produced uncertainty. The pentaquark is produced and recon- 
structed with momentum greater than 25 GeV/c. 

the total likelihood integral to obtain the 95% CL upper limit. To obtain the 
limit on cross section requires a different approach because of the significant 
systematic uncertainties. We use a method based on a note by Convery [24] 
which is inspired by the Cousins and Highland [25] philosophy for including 
systematic uncertainties. The Cousin and Highland prescription is appropri- 
ate for low background experiments with Poisson errors while the Convery 
proposal is applicable to the Gaussian errors which result from the large back- 
ground in our case. Modifications to the Convery approach are made to give 
an exact solution [26] when we include efficiency uncertainty systematics. 

If systematics are not considered, an analysis using a maximum likelihood fit 
returns a central value for the branching ratio (B) and a statistical error (a B ). 
The likelihood function is 



p(B) oc exp 



[B - Bf 
2a% 



(2) 



Including the uncertainty on the efficiency (a e ) changes the likelihood to: 



p(B) oc 



B 2 



1 

7& 



exp 



— (B — B) 2 



erf 



erf 



BBai 



a 



V2a e a B ^jB 2 a 2 

(S -l)a 2 B - Ba 2 e (B - BS) 
V2Sa e a B JB 2 a 2 + a 2 B 



■ (3) 



We integrate Eq. 3 from to oo to obtain the total probability and then 
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integrate from to the point at which 95% of the total probability is included 
to obtain the 95% CL upper limit. The branching ratio B of Eq. 3 is simply 
the relative cross section times the unknown pentaquark branching ratio as in 
Eq. 1. The relative uncertainties on the efficiency for the signal and normalizing 
mode are added in quadrature to become a e in Eq. 3. Furthermore, S is the 
relative efficiency between the signal and normalizing modes and Ob is the 
statistical uncertainty on the branching ratio due simply to the uncertainty in 
the signal yield. 

Figure 12 shows the results for °'^ 5 ^rfaoWT' ~ ' with an assumed natural 
width of (15) MeV/c 2 for the top (bottom) plot. This is the result corrected 
for all undetected particles. The shaded band shows the ±lcx limits with statis- 
tical uncertainties only; the line in the middle of the band is the central value. 
The top curve shows the 95% CL upper limit using the method described 
above including statistical and systematic uncertainties. The curve between 
the full upper limit and the la band is the 95% CL upper limit using the 
method described above with no systematic uncertainties included. The large 
systematic uncertainties are due to the attempt to correct for the significant 
fraction of particles outside of our acceptance. While this systematic uncer- 
tainty significantly degrades the limit, the production times branching ratio of 
the pentaquark relative to S*(1530)° production is still less than 0.032 (0.091) 
at 95% CL over the mass range 1.5 to 2.4 GeV/c 2 for a natural width of 
(15) MeV/c 2 . The background, as shown in Figure 7, is rising so rapidly below 
1.5 GeV/c 2 that fitting becomes very sensitive to the form assumed for the 
background. Consequently we do not quote upper limits below 1.5 GeV/c 2 . 

The plots in Figure 13 show the same results for the restricted range of mo- 
mentum greater than 25 GeV/c. That is, they show limits on relative cross 
sections for particles (Sg" ~, S*(1530)°) produced with p > 25 GeV/c. 



4 Conclusions 



We find no evidence for pentaquarks decaying to S ir in the mass range of 
1480 MeV/c 2 to 2400 MeV/c 2 . In contrast, we observe about 65 000 H*(1530)° - 
2~7r + particles which have a very similar topology and energy release. We set 
95% CL upper limits on the yield over the entire mass range with a maxi- 
mum of 600 (3000) events for an assumed natural width of (15) MeV/c 2 . 
We also obtain 95% CL upper limits on the cross section for pentaquark pro- 
duction times the branching ratio to S~7r~ relative to S*(1530)° — > E~n + . 
These limits are determined for two cases. The first case is for parent par- 
ticles produced at any momenta where we find a maximum upper limit of 

— a'ff(i530)o) B ~ ,r "^ < °- 032 (°- 091 ) at 95 % CL for a natural width of 
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Fig. 12. — ^ — o-( P *(i530)°)" versus mass. Top (bottom) plots show results for 

a natural width of (15 MeV/c 2 ). The shaded region encompasses the la 
statistical uncertainty with the central value in the middle. The top curve shows 
the 95% CL upper limit including systematic uncertainties while the middle curve 
is the 95% CL upper limit with statistical uncertainties only. 

(15) MeV/c 2 . In the second case we measure the relative cross sections for 
parent particles with momenta above 25 GeV/c (a region of good acceptance) 

and calculate 95% CL limits of ^ l^i^o^^ < ' 007 (°- 019 ) for a 
natural width of (15) MeV/c 2 . 

The only experiment reporting an observation of the is NA49 [7] which 
shows about 15 S*(1530)° — > H~7r + candidates, while reconstructing 38 — > 
2~7r~ candidates. The FOCUS results for photon interactions presented here 
represent samples that are more than 4000 times larger, show no evidence for 
a state S5 ~ — > H - 7r~, and are in marked contrast with the NA49 results for 
pp interactions. 
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Fig. 13. — ^ — ct(-*(1530)")" ^ or P > 25 GeV/c versus mass. Top (bottom) plots 

show results for a natural width of (15 MeV/c 2 ). The shaded region en- 
compasses the la statistical uncertainty with the central value in the middle. The 
top curve shows the 95% CL upper limit including systematic uncertainties and is 
virtually indistinguishable from the middle curve which shows the 95% CL upper 
limit with statistical uncertainties only. 
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